Quantum nondemolition-like, fast measurement scheme for a superconducting qubit 
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We present a measurement protocol for a flux qubit coupled to a dc-Superconducting QUantum 
Interference Device (SQUID), representative of any two-state system with a controllable coupling to 
an harmonic oscillator quadrature, which consists of two steps. First, the qubit state is imprinted 
onto the SQUID via a very short and strong interaction. We show that at the end of this step the 
qubit dephases completely, although the perturbation of the measured qubit observable during this 
step is weak. In the second step, information about the qubit is extracted by measuring the SQUID. 
This step can have arbitrarily long duration, since it no longer induces qubit errors. 

PACS numbers: 03.65.Yz, 85.25.Cp, 03.67.Lx, 42.50.Pq 



I. INTRODUCTION 

The quantum measurement postulate is often viewed 
as the most intriguing assumption of quantum physics. 
Much of it has been demystified by the study of the 
physics of quantum measurements. The dynamics of 
the measurement process can be described by a coupled 
many-body Hamiltonian, consisting of the system to be 
measured and the detector with a heat bath component 
[TJ [2J. Thus, the measurement process can be investi- 
gated using the established tools of quantum mechanics 
of open systems [31 SI El E] ■ 

Most interest has been focused on the physics of weak 
measurements, where the system-observer coupling can 
be treated within perturbation theory. Famously, this 
research has shown that only a certain class of measure- 
ments satisfy von Neumann's quantum measurement pos- 
tulate [3 [5] and indeed project the system wavefunction 
onto an eigenstate of the measured observable. Measure- 
ments of this type are termed quantum nondcmolition 
(QND) measurements. Within the weak measurement 
paradigm, the QND regime is achieved when the mea- 
sured observable is a constant of the free motion and 
commutes with the system-detector coupling Hamilto- 
nian. Weak QND measurements have been investigated 
in various systems, ranging from spins to oscillators and 
even photons 0[Tia[TTJ[T21EllIl|TSl[TSlin]- 

The dynamics of the weak measurement process has 
practical relevance in the context of quantum computing. 
Specifically, superconducting qubits have been proposed 
as building blocks of a scalable quantum computer |181 
[19j |20l [21] , and a fast measurement with a high resolution 
and visibility is important for readout and also for error 
correction. 

There are a variety of different measurement tech- 
niques used in superconducting qubits. Weak measure- 
ments can be performed using single-electron transistors 
[18] . A different approach is the switching measurement, 
where the detector switches out of a metastable state de- 



pending on the state of the qubit [22 ] 123 ] 124 ] 125 ] 126]. 
Such switching measurements have been a quite success- 
ful readout scheme for many superconducting qubit ex- 
periments to date. However, the dissipative nature of 
the switching process imposes limitations on the mea- 
surement speed and perturbs the qubit state. 

A QND measurement could be achieved by using a 
pointer system, and measuring one of its observables 
influenced by the state of the qubit [27]. Recent de- 
velopments of such detection schemes, using an oscil- 
lator as the pointer, have led to vast improvements 
[H HO EH Ell E3 E3 EI] over previous measurement 
protocols. 

It has previously been shown [35j ESI E3 EH] that in- 
finitesimally short interaction between a qubit and an 
oscillator is sufficient to imprint information about the 
state of the oscillator onto the qubit. The similar idea 
of using a short interaction to transfer information about 
the qubit into the oscillator has been used [39] m a dis- 
persive readout scheme. In this case, after a short in- 
teraction, the state of the oscillator contains information 
about the qubit which can be extracted by further mea- 
suring one of its observables, for example, momentum. 
However, this scheme did not take possible bit flip errors 
into account. These errors may occur in the short yet 
finite time when the qubit is in contact with its envi- 
ronment. Thus, the full power of a quasi-instantaneous 
measurement has not yet been explored. 

In this paper we describe the effect of an ideally ex- 
tremely short and arbitrarily strong interaction of a qubit 
with its environment (consisting of a weakly damped har- 
monic oscillator). We investigate the back-action on the 
qubit when the measured observable does not commute 
with the Hamiltonian describing the interaction with the 
environment, and study how close this result approxi- 
mates the QND measurement. 

We study a setup consisting of a flux qubit inductively 
coupled to a dc-SQUID magnetometer. The flux qubit 
consists of a superconducting loop with three Josephson 



junctions [301 EJJ- For flux bias near odd half-integer 
multiple of h/2e, the qubit is represented by two circu- 
lating current states with opposite directions. During 
the entire measurement process the SQUID is coupled 
to measurement circuitry, with associated dissipative el- 
ements. However, it never switches out of the zero dc- 
voltage state. The qubit-SQUID interaction of arbitrary 
strength is turned on only for a short time by applying 
a very short bias current pulse to the latter. During this 
time, information about the qubit is imprinted onto the 
SQUID and can later be extracted from it during the 
post-interaction phase by monitoring voltage oscillations 
across the device. When the current pulse is switched off, 
the qubit-SQUID interaction ideally vanishes and the en- 
vironment no longer perturbs the qubit. Thus, one can 
afford a long time to measure the SQUID and determine 
the state of the qubit. 

In section [TTJ following Ref. [22] , we model the qubit- 
SQUID system by a two-level system linearly coupled 
to a dissipative oscillator. We describe the evolution of 
this system by means of a master equation in the Born- 
Markov approximation [43] . valid for the underdamped 



SQUID. In section [ITT] we discuss the qubit-oscillator 
evolution during both interaction and post-interaction 
phases. We study the qubit dephasing and relaxation 
during the interaction phase. We show that, at the end 
of this phase, the qubit appears completely dephased. In 
other words, the qubit has been measured and its infor- 
mation has been transferred in the form of a classical 
probability to the oscillator. During the same time inter- 
val, we find that qubit relaxation has remained negligible. 
For the post-interaction phase we describe the evolution 
of the oscillator under the influence of the environment, 
starting from the state prepared by the interaction with 
the qubit. Technically, extracting the qubit information 
amounts to measuring the amplitude of the ringdown of 
the oscillator momentum. In section IIVI we discuss some 
of the details involved with implementing this measure- 
ment scheme. 



II. MODEL AND METHOD 

We study a flux qubit inductively coupled to a dc- 
SQUID, with one possible setup shown schematically in 
Fig. [l](a). We describe a more detailed setup for imple- 
menting this scheme in section |TV] 

The SQUID is characterized by a two-dimensional 
washboard potential for the two independent phases cor- 
responding to the two junctions [44] . Their sum couples 
to bias current driven through the SQUID, while the dif- 
ference of phases couples to the magnetic flux applied 
to the SQUID. The small oscillations in these two direc- 
tions can have vastly different characteristic frequencies. 
In particular, a small geometric inductance and a low 
critical current can make the flux mode frequency large 
while a shunt capacitor can lower the bias current mode 
frequency substantially. In the limit of very different fre- 




(b) 



FIG. 1: (Color online), (a) Simplified circuit consisting of a 
flux qubit inductively coupled to a SQUID with two identical 
junctions and shunt capacitance C's- The SQUID is driven 
by an bias step-like dc pulse -Zs(t) and the voltage drop V(t) 
is measured by a device with internal resistance R. (b) Il- 
lustration of the measurement scheme: coupling (t — 0) and 
decoupling (t — r) of the qubit and the SQUID (oscillator) 
and the evolution of a point of mass in the transition of po- 
tential from one harmonic oscillator to a superposition of two 
displaced oscillators and back. The dashed (red) and the con- 
tinuous (green) lines correspond to the different states of the 
qubit. 



quencies, one can approximate the SQUID dynamics as 
that of a one-dimensional oscillator in the bias current 
direction, with the position of the oscillator minimum 
dependent on both Ig and the total flux coupled to the 
SQUID which, for example, could vary depending on the 
state of the qubit. 

The setup of Fig. [lj (a) can be described by the effective 
Hamiltonian [32] 



H = Hs + Hj + Hb, 
H s = hwa z + hda x + hn(a^ a + 1/2) 

+ h(e(t) -Q{t-T)){a + a)){a zl + K), 



(1) 



where H$ is the Hamiltonian for the qubit-SQUID oscil- 
lator system, w is the qubit energy and 6 the tunneling 
matrix element, Hb is the Hamiltonian for the dissipative 
environment of the measurement circuitry, Hi describes 
the interaction between the SQUID oscillator and the 
environment, and is the Heaviside step function. We 
note that for a continuous shape of the current pulse sim- 
ilar results are expected, as long as the switching is not 
adiabatic. 

Here the SQUID is described, in the lowest-order ap- 
proximation, by a harmonic oscillator with frequency f2, 
i.e. the plasma frequency of the bias current degree of 
freedom. This frequency also depends on the applied 
bias current, as shown in Appendix [5] This dependence 
leads to a enhanced ring-down frequency after the pulse 
is switched off. This change in the SQUID plasma fre- 
quency does not, in the first approximation, depend on 
the qubit state, therefore it will not qualitatively affect 
this method of discrimination. For the following deriva- 
tion we assume the SQUID plasma frequency constant 
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(the value during the bias pulse), noting that the ring- 
down oscillations occurring in the post-interaction phase 
have in practice a somewhat higher frequency, but oth- 
erwise unchanged behavior. 

The dispersive, next-to-leading order component of the 
qubit-oscillator coupling (39] becomes significant in the 
absence of a linear component for very weak bias pulse, 
which is not the limit we investigate here. In the follow- 
ing, the effects of the linear component are investigated. 
We focus on the regime where the qubit-SQUID inter- 
action displaces the state by more than its zero-point 
fluctuation but does not yet explore the classical nonlin- 
earity. The first consequence of the nonlinear component 
may be to add more phase shift to the ringdown oscil- 
lations. In the measurement protocol proposed here we 
assume a symmetric SQUID. 

The qubit-oscillator coupling strength is tuned by the 
bias current Ib (IS]. When Ib = 0, the qubit and the 
SQUID are decoupled. By using a fast current pulse, 
the qubit-oscillator interaction of arbitrary strength 7 is 
turned on only for the short time r allowing information 
about the qubit to be imprinted onto the oscillator. Dur- 
ing this time, the SQUID oscillator is displaced according 
to the qubit state. After the coupling is switched off, the 
SQUID oscillator phase particle returns to the original 
position after undergoing ring-down oscillations that de- 
cay with a damping determined by the SQUID measure- 
ment circuitry. The parameter K describes the strength 
of the bias current kick induced in the oscillator, caused 
by the abrupt shift in the minimum of the SQUID po- 
tential energy from the bias current pulse, in the absence 
of a qubit. For the expressions of the parameters 7 and 
K and their explicit dependence on Is see Appendix |B| 

During the entire measurement process the oscillator 
is coupled via a linear Hamiltonian Hj 



Hr = 



E 



to a dissipative environment described by a bath of har- 
monic oscillators 



— HuJi I b\ b 



(3) 



with Ohmic spectral density J(lu) — Xfh5{uj — u>i) — 
mhKUj<d(oj — oj c ) /tt [46]. Here [k] — s _1 is the photon loss 
rate. The cut-off frequency lu c is physically motivated by 
the high-frequency filter introduced by the capacitors. 
This environment represents the dissipative element con- 
tained in any measuring device. 

We now describe the dynamics of the qubit and SQUID 
oscillator during the various phases of our measurement 
scheme. 



A. The interaction phase 

At t = 0, before the bias current is rapidly pulsed on 
and the qubit and SQUID interact strongly, we assume 



the factorized initial state p(0) = /3s (0) ®/5b(0). The os- 
cillator interaction with the bath is supposed to be weak, 
and assuming a Markovian environment, we obtain the 
standard master equation for the qubit-oscillator reduced 
density matrix ps(t) — Ttb {p(t)} in the Born approxi- 
mation 



Ps{t) 



1 r 

ih 

1 

W 



H s ,Ps(t) 
t 

dt'TiB 





(4) 



Hr 



Hi(t,i?),p s (t)®p B {0) 




where 

Hj{t,t') = , Ut = Texp 



(5) 

and T is the time-ordering operator. 

This approach is valid at finite temperatures /;bT > 
Hn and times f > 1 /uj c [3 , 47 , which is the limit we will 
discuss henceforth. 

In the qubit a z eigen-basis the density matrix and the 
qubit-oscillator Hamiltonian read 



Ps 

Hsaa 



(G) 



ptt Pn 
v Ait Au 

= H sn - HS, r a = (a\*,\a), a G {U}, (7) 
= h(r a w + VL{a i a+ 1/2) (8) 
+ (rvy + 1*0(0+0*)). 



In the following, we assume that the environment acts 
on each matrix element of ^ in the same way. This is a 
valid assumption in the case of very weak damping and 
8/w <C 1 for an Ohmic bath. Within this assumption we 
obtain 



Paa - 

hi : 
pit : 

where 

£P<T<T> 



1 

ih 
1 

ih 
1 

ih 



[H a<J ,p aa ] - iSr a (p l1 - p u ) + £p aa , (9) 
(H n p n - p n H n ) + iS(p n - p u ) + Cp n 

( H iiPn - Pn H ii) - iS (Pn - hi) + £ht> 



-K(a)ap aa > 
2nn{a> ap aa 



+ p aa >a)a - 2ap (TCT /a t ) 



'da' — ap 



(10) 
a) p aa ia). 



At t = we assume a factorized initial state for the 
qubit-oscillator reduced density matrix 



Ps(0) =/5 q (0)®/5 H o(0), 



(11) 



and use the Wigner representation of the oscillator den- 
sity matrix in phase-space |48j 



1 



/5ho(0) = - / daXo(a)D(-a), 
D(—a) = exp {-~aa) + a*a) , 



(12) 
(13) 
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where xo is the Fourier transform of the Wigner function. 
We assume the oscillator to be initially in a thermal state 



Xo 



(a) = i-exp(-||a| 2 ), r, = l + 2n(Q), (14) 



where n(f2) is the Bose function at bath temperature T. 
The qubit is assumed to be initially in the pure state 
I*) = (ftl t) + 9i<e i0 l I) such that 



pq(0) = U*> 



(15) 



For the corresponding Wigner characteristic functions 
we obtain the following coupled partial differential equa- 
tions: 

Xaa = (i(r a j + K)(a + a*) +in(ad a - a*d a *) 

+ fk^-v^xit-xu)). ( 16 ) 

Xti = (2i7(9 a » — 9 a ) + iQ(ad a — a*d a *) — 2iw 

+ iK(a + a*) + £>)xu -i<Kxu -Xn))> 
Xj.t = (-2i7(9 Q . - 3 Q ) + iri(a(9 a - a*9 Q .) + 2iw; 
+ iif(a + a!*) + 2?)X1T + i( KXU - Xtt)), 
where the differential operator T> is given by 

V = -K(ad a + a*d a .) -r]K\a\ 2 . (17) 

To solve these equations, we approximate the inhomogc- 
neous parts, in the limit of short time r and weak tun- 
neling S, by 

Xro :(()-v(0) + iv(0), a,(/e{t,i}. (18) 
For details on the solution see Appendix [X] 

B. The post-interaction phase 

The state prepared by the interaction with the qubit 
at t = t, as the bias current pulse ends, is described by 

P(r) = W)W™>(t)®Pb(Q)- (19) 

<T,a'e{u} 

Since the system Hamiltonian no longer contains any 
qubit-oscillator interaction, we can write the time evolu- 
tion of this density matrix as follows 

p{t) = U q (t)\a){a'\Ul(t) 
<T,ff'e{U} 

(O)^o-bW. ( 20 ) 
where the evolution operators are given by 

U q {t) = eM-i{t-r)(Sv x + wa z )), (21) 
Uho-b = Texp I / dt' — '- .(22) 



In the reduced density matrix 
p s {t) = Tv BP {t)= Y, U q {t)\a){a'\fj\{t) (23) 



Tr, 



{ UnO-B (t)p aa , (t) <g> p B (0) l>HO-B CO } 



we can treat the time evolution of the oscillator compo- 
nents by means of a master equation in the Born-Markov 
approximation and, in a similar manner to Eq. @, we 
obtain 



Paa-'(t) = -i£l[a*a,i 



(t)} (24) 
Hi,[Hi(t,t'),p aa ,(t) ® p B (0)] 



Using the Wigner representation 
1 



Paa> (t) 



d a Xcra'{ot,t)D{-a), 



(25) 



we obtain the differential equation 

Xaa'(a,t) = (ifl(ad a ~ a*d a .) + V)xaa'(a,t),(26) 

with the initial condition prepared at the end of the in- 
teraction phase 



Xaa'(a,r) = Xaa'(a,T), 

and the analytic solution 

W(«,t) = w(ae" (t " r)(,C ""V) 



x exp 



(|i«i 2 (. 



-2(t-r)« 



- 1 



(27) 



(28) 



The reduced density matrix in the post-interaction phase 
is given by 



hit) 



d 2 a Xs s>(a,t)D(-a), (29) 



where 



Xss>(a,t) =Y (s\U q (t)\a)(a'\Ul(t)\ S ')x**>(a,t).(30) 
^'e{T,|} 

In the the post-interaction phase, the qubit and the oscil- 
lator are decoupled. The trace of the oscillator-bath part 



in Eq. (20) is time independent, as one can see after a 



circular permutation of the involved operators. One finds 
that the qubit time evolution is given only by the unitary 
U q , and thus is independent of the oscillator. Physically, 
this means that in the post-interaction phase no further 
information about the qubit can be transferred to the 
oscillator-bath system, and thus the qubit suffers no fur- 
ther decoherence. 



III. RESULTS 

In this section we analyze the qubit decoherence and 
the evolution of its detector, the dissipative oscillator, 
during the entire measurement process. 
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A. Qubit decoherence 

During the interaction phase, t £ (0, r), the qubit is in 
contact with an environment represented by the dissipa- 
tive oscillator, and thus subject to decoherence. 

The qubit can be prepared in a well defined state by 
thermal relaxation or (if the temperature is too high) by 
measurement post-selection and conditional rotation by 
microwave pulses. 

We analyze the qubit relaxation described by 

(*.)(*) =47r(Xn(0>*)-Xu(0.*)) ( 31 ) 
and from Eq. ( |A5[ ) we obtain the analytic result 

(*,)(*) = (q^-ql)(l-2t 2 6 2 ) (32) 
+ 4qiqit5(twcos(4>) +sm((j))). 

We observe that the above expression is identical with 
the expansion up to the second order in time of {a z ){t) 
when the qubit evolves under the free Hamiltonian H q 
only. Thus, the evolution of (a z )(t) in this short time 
expansion is indistinguishable from the free evolution of 
the unperturbed qubit. This can be understood as fol- 
lows: the observable a z commutes with the environment 
coupling, but is not an integral of the free motion, as 
required for a QND measurement [T]. Thus, the pertur- 
bation of the measured observable comes only from the 
free evolution of the system. One can restrict this per- 
turbation by reducing the time r when it takes place. 
Fig. [2] (a) shows the evolution of (a z )(t) for a set of pa- 
rameters closely related to a feasible experiment, see also 
Appendix [b] The initial state chosen for panel (a) was 

Furthermore, we analyze the qubit coherence (a x ) 
which is given by 



(a x ){t) = 87rRex n (0,t), 



(33) 



and can be evaluated from Eqs. ( A13|A14 ), where 
X 1 |" h (0,t) can be integrated numerically. 

We observe that, if the interaction time r is long 
enough to allow the oscillator a full period evolution, one 
finds a revival in the qubit coherence at the end of this 
period. As the oscillator returns to (almost) its initial 
state, the information about the qubit is "erased" from 
the oscillator, as the oscillator states corresponding to 
I |) and I J.) arc no longer discernible. The height of 
the coherence revival peaks at fit = 2im decays in time 
as the information about the coupled qubit-oscillator sys- 
tem flows (irreversibly in this case) into the environment. 

The qubit dephasing for the same parameters of Ap- 
pendix [b] is shown in Fig. [2](b). The appropriate initial 
state for this study is the equal superposition (l/v / 2)(| t 
) + I [))■ We observe that the qubit appears completely 
dephased after the strong interaction with the damped 
oscillator, such that only a classical probability is im- 
printed onto the latter. 




0.4 0.6 

t/x 



FIG. 2: (a) Evolution of (a z ) with qubit initially in | f) state, 
(b) Dephasing from the 1/V2(| T) + 1}) state for the time 
r that the qubit is in contact with the oscillator. For both 
plots, the following parameters were used: Q/(2n) — 0.97 
GHz, k/O = 10" 2 , w = O, Q.T = 1.83, St = 0.015, jt = 3, 
T = 30 mK. The assumed values of the circuit parameters are 
given in Appendix |B| 



In Fig. [2] (a) we observe that the relaxation from the 
excited qubit state is very weak during the interaction 
time, as (<x z ) differs at most by 10 -3 from the initial value 
of 1. This combination of low coherence (b), indicating 
the fact that the information about the qubit has been 
imprinted onto the oscillator, and very low relaxation (a) 
demonstrates that the first step of the measurement pro- 
tocol produces a good starting point for the second one, 
the oscillator readout. The negligible relaxation brings 
the qubit close to QND dynamics. 

We observe that the qubit coherence time is essen- 
tially dominated by the coupling between the qubit and 
its complex environment 7 -1 such that it is desirable to 
achieve 77 > 1. The relaxation of the qubit has been 
described in the first order in time, and essential to the 
almost-QND result is that tS -C 1. We note that the 
implied condition 7 <C S contradicts none of our approx- 
imations, and can also be realized in experiment. 



B. Detector dynamics 

In this section we study the evolution of the damped 
oscillator, which represents the detector. To achieve the 
strong qubit-oscillator coupling during the short interac- 
tion phase required to imprint the qubit state onto the 
oscillator, one needs a bias current pulse that approaches 
the critical current for the SQUID. Nonetheless, it is im- 
portant that the SQUID does not switch out to the run- 
ning state during the bias current pulse. For the parame- 
ters given in Appendix[B] we can evaluate the SQUID es- 
cape rate [39] from the zero- voltage state during the bias 
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current pulse in the regime of quantum assisted thermal 
activation (fc^T > Ml) 



sinh 



r — 

- 1 - SW 



\2k B T ) 



exp 



-At/ 



(34) 



where AU is the potential barrier. We obtain, for the 
worst case, T sw « 3.6 • 10 7 s _1 such that the escape time 
is much larger than the duration of the bias current pulse. 

The output of the detector is the time dependent volt- 
age across the SQUID, which is proportional to the mo- 
mentum of the oscillator. The probability distribution of 
momentum is given by 



P(p,T,t) = n(5(p-p)) 



(35) 



2 / da x ^ Xcrcr(a x ,t)exp 



T6{U} 



V = 



mflh 



- tot 



(36) 



where, in the post-interaction phase (t > r), X(ra( a xit) 
also depends on r via its initial condition. The expecta- 
tion values for the n th moment of the oscillator momen- 
tum and position are then 



(P n )(t) 



(x n )(t) 



(-l)"aj n V^K,t)| a ^o,(37) 



/ j A.tct\ 

<re{T,l} 



V / <re{U} 



Furthermore, in the post-interaction phase we have, from 
Eq. pi 



E 



Xaa = ]T (<j\U q (t)\ S )( S '\Ul(t)\a) X 
s.s'.aefT.i} 



E^ ss ' 



(38) 



which shows, as expected, that no measurement of 
the oscillator can provide information about the post- 
interaction evolution of the qubit, provided this evolu- 
tion is unitary (i.e. the qubit is not being measured by 
something else). 



For the evaluation o f bo th Eqs. (35 37 1 the s 



integration in x™„, Eq. (A7l should be evaluated last 



Thus, one obtains an analytic (but rather long) expres- 
sion for the expectation value of momentum, while for the 
probability density a numerical s-integration is required. 
Nevertheless, the components originating in x^J 11 turn 
out to be dominant, and we give their analytic expres- 
sions in the following: 



(p)(r,t) = (p)hom(r,t) + (p)inh(r,t) 

(p)hom 



(39) 



(K + 7 <Zt-7<7i) ^~ {t - T)K 



-(t--r)in 



—k — iSl 



+ <E (t-r)iO; 



-r(K-iO) 



The explicit form of the probability distribution of mo- 
mentum, Eq.d35|, is given by 



where 



P(p,T,t) = P hom (p,T,t) + P inh ( P ,T,t), (40) 



iW P , r, t) = V ^ exp - 



-(t-T)if2 



to 1 - ( E - T (' t + if2 ) 



+ <E (t-r)tn : 
k — in —k + iCl 



The results above refer to the post-interaction phase 
t > t. For the interaction phase, t £ (0,r), the prob- 
ability distribution of momentum is given by P(p, t, t) 
in Eq. (40 1 and the expectation value of momentum by 
(p)(t,t) in Eq. (39), i.e. by replacing r by t. 



The expectation value of momentum (p)(r, t) in the 
post-interaction phase contains information about the 
qubit initial state. We observe that the momentum os- 
cillations corresponding to the two different initial qubit 
states | t) and | |) for t > r are in phase. Disregarding 
the inhomogeneous contributions, which are relatively 
small in the limit of small tS, the envelope of the ho- 
mogeneous part is given by 



Vk 2 + ft 2 



• v/-2ce- Kr cos(rfi) + <sr 2KT + 1. (41) 

Fig. |3] illustrates the phase-space trajectories of the os- 
cillator corresponding to the qubit being in either the 
| t) ° r I I) state. During the interaction phase the sys- 
tem moves away from the origin. After switching off the 
interaction, the trajectories spiral back towards the ori- 
gin, without crossing. For K = the trajectories are 
symmetric with respect to the origin, while intro- 
duces an asymmetry. We note that the artificial situation 
K = includes only the bare oscillator response for the 
different qubit states. This situation has been introduced 
in order to more easily illustrate the difference between 
the two oscillations. 

Fig. [4] shows the output of the detector for the two 
qubit states | J.) and | j). 

The standard condition for the possibility of single- 
shot readout, i. e., the maximal separation of the two 
peaks corresponding to different qubit states in the prob- 
ability distribution Eq. ( 35 1 should be larger than the 
peak width, is given by 



|.4(1,0) -.4(0,1)1 



> 1, 



(42) 



where the envelope (41 1 has been evaluated at t — t. We 



note that and q[ are continuous variables with values 
between and 1 and the condition presented above takes 
into account the extremal case of the difference between 
the states | |) and | J,). The result is independent of K. 
For the parameters of Fig. [4] we have e ~ 2.5. 
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FIG. 3: (Color online). Phase space representation of the 
oscillator trajectories ((x)(t), (p)(t),t) corresponding to the 
two qubit states | J.) (dashed, red) and | |) (continuous, green) 
for the parameters given in Appendix [B] an oscillator quality 
factor of 10, with K = (a) and K 7^ (b). Projections on the 
(x,p), (x,t) and (p, t) planes are included. Both trajectories 
start at the origin and move away from it under the influence 
of the interaction with the qubit. At the point marked with • 
the interaction is switched off, and the system evolves freely 
spiraling around the origin. The trajectories circle around 
each other without crossing. 



the same phase and amplitude as those of the measure- 
ment SQUID oscillator for one of the two qubit states. In 
this case, the resultant signal after the subtraction would 
be exactly zero for perfect cancellation when the qubit 
state causes the two SQUID oscillators to have identical 
ringdown signals. A residual ringdown oscillation would 
be produced for the other qubit state. This scheme re- 
quires that the two SQUIDs receive an identical kick and 
begin their ringdown oscillations at the same time. This 
can be achieved by splitting the bias current pulse sig- 
nal along two separate lines, one going to each SQUID, 
as shown in Fig. [5] where the layout is such that the 
reference SQUID has a vanishing coupling to the qubit. 



SQUID oscillator Gradiometer HEMT 

& qubit : Microstrip SQUID 

1 1 1 1 amplifier 




reference SQUID i3Q_mK , 



oscillator 



O: 2 

CD 

cd 1 

-4— * 1 
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FIG. 4: (Color online) Probability distribution of output volt- 
age (density plot, dark color indicates high and white low 
density) and expectation value of momentum for the two 
qubit states | J.) (dashed, red) and | f) (continuous, green). 
Here fi/(2vr) = 0.97 GHz, 0,/k = 20, w = ft, fir = 1.83, 
St = 0.015, 7T = 3, T = 30 mK. The assumed values of the 
circuit parameters are given in Appendix [B] 



IV. PRACTICAL IMPLEMENTATION 

A possible measurement protocol involves discriminat- 
ing the amplitudes of the ringdown oscillations corre- 
sponding to different qubit states. As demonstrated by 
Eq. (41 1, the amplitude difference is independent of K. 



This discrimination could be performed more accurately 
with an interferometric technique, where ringdown oscil- 
lations from a second, reference SQUID oscillator that is 
not coupled to the qubit are combined with those from 
the original SQUID oscillator. The reference SQUID is bi- 
ased such that it undergoes ringdown oscillations with 



FIG. 5: Circuit diagram for SQUID oscillator and qubit, 
along with reference SQUID oscillator, dual-input gradiome- 
ter microstrip amplifier and a cryogenic High Electron Mobil- 
ity Transistor (HEMT). Dashed boxes indicate different chips 
and/or different temperatures. 

Fig. [6] shows the total signal, i.e. the difference of the 
ringdown oscillations from the measurement and refer- 
ence SQUIDs for the two qubit states. We have consid- 
ered the case where the total flux bias for the reference 
SQUID is equal to the total flux bias for the measure- 
ment SQUID in the case where the qubit state is | j). 
In this case the difference signal is smeared around for 
the qubit in state | |). If the qubit is in the | J.) state, 
the output signal oscillates with an amplitude is given 
by the difference between the two ringdown oscillations 
in Fig. [4] 

The subtraction of the two ringdown signals can be 
achieved by using a microstrip SQUID amplifier arranged 
as a gradiometer with two separate microstrip inputs 
with their senses indicated in Fig. [5] [50] . The microstrip 
SQUID amplifier consists of a dc SQUID with a multi- 
turn superconducting input coil above a conventional 
SQUID washer, where the signal is connected between 
one side of the input coil and ground and the other end of 
the input coil is left open. Input signals near the striplinc 
resonance frequency, related to the total length of the in- 
put coil, typically of the order of 1 GHz, couple strongly 
to the SQUID loop and the SQUID produces an output 
signal with a gain of ~ 10 — 20 dB [5T]. A gradiometer 
microstrip SQUID amplifier for amplifying the difference 
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FIG. 6: (Color online) Probability distribution of output volt- 
age (density plot, dark color indicates high and white low den- 
sity) and expectation value of momentum for the two qubit 
states | I) (dashed, red) and | T) (continuous, green). Here the 
contribution of the reference SQUID has been introduced. Pa- 
rameters: Cl/(2ir) = 0.97 GHz, Q/k = 20, w = Q, Qt = 1.83, 
St = 0.015, 7r = 3, T = 30 mK. The assumed values of the 
circuit parameters are given in Appendix [B| 



between two separate signals near the stripline resonance 
can be produced as a straightforward extension from pre- 
vious microstrip SQUID layouts by using a SQUID geom- 
etry with two loops and a separate stripline coil coupled 
to each of the loops, with one signal input connected to 
each stripline [SU]. With no crosstalk between the two 
inputs, the circulating currents in the two loops of the 
SQUID amplifier cancel out when the input signals are 
identical, resulting in a vanishing output signal. Thus, 
with the arrangement in Fig. [5j the microstrip SQUID 
amplifier produces the difference between the two oscil- 
lator ringdowns. Of course, in any practical gradiometer, 
there will be non-zero crosstalk, where a signal at one in- 
put induces circulating currents in the other loop of the 
SQUID amplifier. However, for reasonable layouts of the 
device, this crosstalk could be kept at the 1% level, thus 
setting a limit on the fidelity of the subtraction [5U] . 

Based on the calculated difference signals for the ring- 
down oscillations in the two qubit states from Fig. [6j one 
must be able to discriminate the oscillations for the | |) 
qubit state from the non-oscillatory signal for the | f) 
state. Thus one needs to resolve a ~ 1 GHz signal with 
an amplitude of ~ 0.5 /iV in a ~ 100 MHz bandwidth, 
i.e. before the ringdown is completed. Microstrip SQUID 
amplifiers operated at 20 mK have achieved noise tem- 
peratures as low as ~ 50 mK [52] . If we assume a con- 
servative noise temperature estimate of 200 mK for our 
gradiometer microstrip SQUID amplifier, this would cor- 
respond to a noise of 250 nV in the 100 MHz bandwidth 
referred back to the SQUID oscillators. Thus, it should 
be possible to discriminate between the two possible out- 
put signals corresponding to the two qubit states in a 
single shot. 

In the non-ideal case, the noise of the reference SQUID 
increases the broadening of the curves in Fig. [6] such that 
the single shot condition ( 42 1 must accommodate another 



width r\. Still, at the parameters used in Fig. [6j this 
condition will still hold. 



V. CONCLUSION 

We have demonstrated that a non-QND Hamiltonian 
can induce a close to QND backaction on the qubit, 
despite arbitrarily strong interaction with the environ- 
ment, provided that the interaction time is very short, 
i.e. the measurement is quasi-instantaneous. The relax- 
ation of the qubit has been described in the first order in 
time and, essential to the almost-QND results presented 
above, is that tS -C 1. 

We observe that the measurement time, i.e., the time 
needed to reduce the qubit density matrix to a classical 
mixture is essentially dominated by the coupling between 
the qubit and its complex environment 7" 1 such that it 
is desirable to achieve jt 3> 1. 

The readout time for the oscillator is restricted only 
by the ring-down of the two possible oscillations of mo- 
mentum, i.e. The amplitude of these oscillations 
is proportional to 7, which again stresses the usefulness 
of a strong qubit-oscillator coupling. If the two peaks 
in P(p, r, t) become separated by significantly more than 
their widths, single shot measurement may become pos- 
sible. 

The method presented above has the advantage of 
a very short interaction between the qubit and its en- 
vironment, compared to e.g. the dispersive readout of 
Ref. [25], and results in a QND-type of readout, with- 
out the requirement of strong, continuous AC driving of 
e.g. Ref. |29j which may induce spurious qubit relaxation. 

As a figure of merit we consider the QND fidelity in 
Ref. [53]. For the parameters used in Fig. [2] and an initial 
qubit state \& = l/\/2(| T) + | |)), our scheme achieves 
at the end of the post interaction phase a QND fidelity 
of 99.92%. 

Furthermore, if the aim is to apply the idea of a short 
interaction with an intermediate system, dispersive mea- 
surement, with all its potential advantages, may be diffi- 
cult due to the continuous driving which implies contin- 
uous interaction between the qubit and its environment. 
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APPENDIX A: SOLUTION FOR THE WIGNER 
CHARACTERISTIC FUNCTIONS 

In this section we solve Eqs. ( 16 ) using the approxima- 
tion ([181. 
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a. The diagonal density matrix elements 



The variable transformation in this case originates from 



We solve the diagonal equations needed for evaluation 
of expectation values such as (p)(t), which characterize 
the output of the detector: 

Xaa = (i( r <r7 + K)(a + a*) + ifl(ad a - a*d a .) 

+ V)x a a-r (y iSxo(a)F(a,t), (Al) 

where 

F(a,t) = 2g T ^sin(^)(i - K(a + a*)t) - 2i(g| - qf)St 

— 2iq^q±cos(<fi)t(r]j(a* — a) — 2w). (A2) 

We perform a variable transformation in order to remove 



the first order derivatives in Eq. ( Al I 



a = 
and obtain 



ZIB s(K-m) , 



= z *<B s( - K+m \ t=s, (A3) 



d sX aa = (i(rV7 + K)<B SK (z<B- sin + z*<s sin ) 

- ? 7K |Z| 2 (E 2SK )X (TCT 

- r !T idxo(z<B s(K - in) )F(z<B s( - K - in \s), (A4) 

which can be solved analytically, and transformed back 
to the initial variables a,t. The solution reads 

XCTCT (a, t) = M^^a, t) - ^X^(«, *), (A5) 



where 



,hom 



(a, t) = cxp 



|a| 2 ?7 



+ i(r CT 7 + K) 



(A6) 



afl-e-*^ 1 ")) o;*(l-(E-'( K + in )) 



■ifi 



■ m 



and 



X ^ h (a, t) = jdsxlT{<*, s)F (a^- in \t - s) 



(A7) 



b. The off-diagonal density matrix elements 



d s a = (— if2 + k)q + 217, 
9 s a* = (ifi + — 217, 



and reads 



a 



k — ifi V / 
2i7 
k — ifi 



f = s. 
We obtain 



(A10) 



9 s xtj. = (— 2iw — r)Ka(z, s)a*(z*, s) 

+ iK(a(z,s) +a*(z*,s)))xn 

- idxo(a(z,s))G(a(z,s),s), (All) 

which can be solved analytically and transformed back 
to a,t. The solution reads 



xn (a,t) = M^ x ft m (M - ^xu («>*), ( A12 ) 



where 



-77 — 2ii«; — 



4^7 (777 k — ilffi) 



4 7 (77;(k 2 - ft 2 ) - 2iX K fi) 

(K 2 + n 2 ) 2 



(A13) 



k + ifi 
K — ji] 
k — ifi 



i(l - e-*( re+in ))a* - 

i(l_ (E -*(«-i")) Q; + 



k + if2 

2(E -t(K-io) 7 

k — ifi 



and 



(A14) 



c f c -(^)q + 2 ( i - c "' ( "" in ^, t -; 

\ UK + fi 



The method and approximations of the previous sec- 
tion can be used to solve the off-diagonal equations. From 
this solution we intend to extract information about the 
qubit coherence (a x )(t). We start with 

Xti = (2i7(9 Q » — d a ) + ifl(ad a — a*d a *) — 2iw 
+ iK(a + a*) + V)xaa-iS X o(a)G(a,t), (A8) 



where 



G(a,t) = qf-qi-tiij-Kiql-qtyia + a*) 

- 4i5g m sin(<£). (A9) 



From the density matrix calculated above we can ex- 
tract information about the qubit relaxation and dephas- 
ing during the short interaction with the dissipative os- 
cillator. 



APPENDIX B: CONVERSION TO CIRCUIT 
PARAMETERS 

In the following we give a recipe [32] to obtain the 
parameters entering the calculation of this paper from 
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the circuit components 



0= vlSC ['-(i? 



2\ 4 



/« =1^1 C S , 



MqsLlB tan (jP m 
7= — — . k 



2tt 
1 



4/x 



tan i 



K 



2RC S ' 

i B r h 



2e V 2mtt 



where $o = h/2e is the magnetic flux quantum for a 
superconductor, Mqg is the qubit-SQUID mutual induc- 
tance, 1% is the effective critical current of the SQUID at 
the particular flux bias, Ib is the amplitude of the dc bias 
pulse applied to the SQUID, Cs the SQUID shunt ca- 
pacitance, R the internal resistance of the measurement 
circuitry, and I q is the circulating current of the local- 
ized states of the qubit. The expression for K is derived 



in the limit of a small geometric inductance, low critical 
current and large shunt capacitor where one can approx- 
imate the SQUID dynamics as that of a single Josephson 
junction with a variable critical current. The momentum 
of the oscillator p and the voltage across the SQUID are 
related by 



V = 



ep 



Csh 1 



(Bl) 



where e is the electron charge. The parameters used to 
generate Figs. |2j [3] gj [6] are 



7 c off = 0.5 • 10~ 6 A, 
C s = 2 • 10~ n F, 
I q = 438 • 10~ 9 A, 



I B = 0.87/f, 

M qS = 100 • 10" 12 H, 

t = 0.3 • 10~ 9 s, 



S/(2n) = 0.8 • 10 7 Hz. 
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